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The calculation of expected spectral line strengths and profiles is a powerful tool for the analysis of the solar 
atmosphere, and other stellar atmospheres. We present here a recipe in seven easy steps for the development of 
such spectral synthesis software. 
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1. Introduction 

The solar spectrum can be calculated by (nu- 
merically) solving the radiative transfer equation: 

Mj- = h ~ Sx (1) 

If we know the source function S\, then finding 
the emergent spectrum is straightforward. With a 
model solar atmosphere, we have the LTE source 
function, S\(z) — B\(z), the Planck radiation 
function, as a function of physical depth. Un- 
fortunately, it is not sufficient to know 5a as a 
function of physical depth z; we must convert the 
physical depth scale to an optical depth scale: 



t a(z) 



p(z r ) Kx (z')dz' 



This requires calculation of the opacity k\; this 
is the core task of a spectral synthesis program. 
The opacity can be considered to consist of a con- 
tinuous opacity, and opacity due to atomic line 
transitions. The main contribution to the contin- 
uous opacity is due to bound-free and free- free H~ 
transitions; other major sources (H^~ molecule, 
photoionisation of metals) also need to be con- 
sidered. The line opacity can be found from the 
line strength and line profile for each line in the 
spectral region of interest. 

2. Model atmospheres 

The solar atmosphere is highly stratified; its 
physical properties vary strongly with height 



1. Calculate S\(z) 

2. Calculate continuous opacity 

3. Calculate line opacities. 
This requires: 

• line strengths 

• damping 

• thermal and mass motion 

4. Calculate optical depth scale 

5. Effect of intermediate scale mass motion 

6. Numerical integration of radiative transfer 
equation JIJ 

7. Non-plane parallel corrections — large scale 
mass motion 

Figure 1. Recipe for spectral synthesis 



within the photosphere. Thus, the photosphere 
can be well represented by a plane-parallel model 
atmosphere. Many model atmospheres are avail- 
able, with those by Holweger and Muller [[] and 
Kurucz [2] being commonly used. 

A model atmosphere tabulates temperature, 
pressure, electron pressure, density, and opacity 
at 5000 A. 

With this data, and assuming LTE conditions, 
then the source function is simply equal to the 
Planck radiation function: 
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3. The continuous opacity 

The continuous opacity (maximum at 8200 A) 
is dominated by the H ion, due to its large inter- 
action cross-section. We can also consider other 
abundant absorbers. If we are only considering 
a small spectral region, we can assume that the 
continuous opacity is uniform. In general, we use 
approximation formulae and table lookup 3 . 

The H~ opacity is due to both bound-free 
(H~ + photon H + e~) and free- free (H + 
e~ + photon <-> H + e - ) processes. Approxima- 
tion formulae for both processes are given by Gin- 
gerich ^ 

Due to its high abundance, we should also con- 
sider bound-free and free-free absorption due to 
neutral hydrogen, which depend on the Gaunt 
factors gu and gm. Approximation formulae for 
gn and gm are given by Mihalas 

Other sources of opacity which can be consid- 
ered include the molecule, for which cross- 
sections have been calculated; a convenient tab- 
ulation for interpolation is available in ,|[. The 
abundances of heavier elements much lower than 
that of hydrogen, and they contribute far less to 
the opacity. Only the most abundant species need 
be considered, Mg and Si being the most impor- 
tant. Calculations by Peach |H] can be used for 
interpolation. 

4. Line opacities 

The opacity due to a spectral line depends 
strongly on the wavelength — which is why it's a 
line — which makes it necessary to calculate the 
line opacity on a sufficiently dense grid of wave- 
length points. For a single (not too strong) line, a 
total wavelength range of « 1 A is enough; a few 
hundred wavelength points will be sufficient. The 
use of a non-uniform grid of wavelength points 
can greatly speed up line opacity calculation, at 
the cost of a greater effort to implement. 

The opacity of a single line depends on 

1. population of the absorber, 

2. intrinsic strength of the transition, 

3. and the line profile. 



The populations of absorbers depend on the 
atomic abundances and the equilibrium between 
the atomic species, the rest of the atmosphere, 
and the radiation field. We can simplify this 
greatly by assuming LTE. 

The strength of the transition — the os- 
cillator strength, or /-value can be ob- 
tained from a suitable source of data such 
as the NIST Atomic Spectra Database at 
http : //physics . nist . gov/ cgi-bin/AtData/main_asd 
or the Vienna Atomic Line Database (VALD) at 
http : //www . astro . univie . ac . at/~vald/ 

4.1. The line profile 

The line profile function of a stationary atom 
is the Lorentz profile: 
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(4) 



Thermal motions will produce a profile that is 
the convolution of the Lorentz profile and the 
Maxwellian Doppler shift profile, giving the Voigt 
profile: 
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where a — AAl/2AA,d and v — (A — Ao)/AA_d, 
where AAi, is the natural (Lorentzian) line width, 
and AAd is the Doppler shift due to the most 
probable speed. 

We can expect a spectral synthesis code to 
spend a major part of the execution time re- 
peatedly calculating the Voigt function. Various 
methods are available |7I8| . Fast methods usually 
use asymptotic methods where possible (since the 
Voigt profile approaches the Lorentz and Maxwell 
profiles in the limit), and other approximate (eg 
power series) methods elsewhere. Alternately, the 
Voigt profile can be calculated by convolution of 
the Lorentz and Maxwell Doppler profiles, using 
Fourier transformation methods for speed. 

But, first, we must find the Lorentzian - the 
natural line width and damping - and Maxwellian 
widths - thermal and other motion. 
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4.2. Damping 

The most important contribution to damping 
is collisions with neutral hydrogen (almost all of 
which is in the ground state). Damping due to 
electrons is 50 times smaller, or even less, and 
damping due to helium about 30 times smaller. 
Other sources can be neglected. 

The line broadening theory developed by 
Anstee and O'Mara |3 has been shown to be ac- 
curate. Collisions with hydrogen atoms in the 
photosphere are fast; the impact approximation 
can be used. The line width (HWHM) is 

AAcou = N / vf(v)a(v)dv (6) 
Jo 

where a is the line broadening cross section, and 
f(v) is the Maxwellian velocity distribution. 

Data calculated using this theory is avail- 
able 513) an d code to calculate broaden- 
ing has been developed (available at 
http://www.astro.uu.se/~barklem) which can 
be used either for stand-alone calculation, or can 
be incorporated into a spectral synthesis pro- 
gram. This code calculates the linewidth per unit 
hydrogen atom density for a given temperature. 

The total Lorentzian width is the sum of the 
collisional damping width, the natural line width, 
and the stimulated absorption/emission width. 
Only the collisional damping width will be im- 
portant in the photosphere. 

4.3. Small-scale mass motion 

We can expect small-scale turbulent motion 
to have a Gaussian velocity distribution (Kol- 
mogorov turbulence) . The effect will be the same 
as that of thermal motion, and the two can be 
combined. The average photospheric microtur- 
bulence is 0.845 km/s. We can simply use this 
average value as a uniform microturbulence, or a 
depth-dependent value can be used. 

4.4. Intermediate-scale mass motion 

Larger scale mass motion, that cannot be 
treated as Gaussian turbulence, will occur. A full 
treatment would require knowledge of the fluid 
flow on length scales of ~ 10 - 100km, and ex- 
tremely difficult problem. In general, such motion 
will act so as to Doppler shift the opacity; this will 
produce broadening if it varies with depth. 



5. Integrating the radiative transfer equa- 
tion 

All of the processes that affect the opacity have 
been considered. The radiative transfer equa- 
tion Q can be integrated. The radiative transfer 
equation Q is well behaved, provided we solve 
it in a reasonable manner. In order to avoid the 
growth of errors, we integrate 

Ix(t = 0)= [ Sx(Tx)eM-T\/»)~dT X . (7) 
Jo M 

In practice, we step through the layers of our 
model atmosphere, from the top to the bottom, 
with the step size for the numerical integration 
being determined by the spacing in the model at- 
mosphere. By the time the bottom of the model 
atmosphere is reached, the contributions to the 
emergent intensity are negligible. 

6. Large-scale mass motion — beyond the 
plane-parallel approximation 

Large-scale flows prevent the photosphere from 
being strictly plane-parallel. Further effects must 
be considered as corrections to the simple plane- 
parallel model. The simplest method is to assume 
that such motions have a Gaussian distribution, 
and will only act to broaden the emergent spec- 
trum (as opposed to the microturbulence, which 
will broaden the line profile at all depths). This, 
however, is an unsatisfactory solution. The main 
motion at this scale is that due to the granula- 
tion. A granular cell is typically about 1000km 
across, consisting of a hot, slowly rising, central 
region, and a cooler, more rapidly falling outer 
region. 

Even a simple model of granulation can greatly 
improve the results of spectral synthesis . If we 
divide the granule into three regions: the rising 
core, the falling outer region, and a transition re- 
gion between the two, we can treat each region 
as purely plane-parallel. We can assume that the 
flow velocity in each region varies with height h 
as 

V = V Q exp(-h/V s ), (8) 

where V s is the velocity scale height. We can 
expect the microturbulence, since it is driven by 
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the large scale flow, to have the same depth de- 
pendence. A macroturbulent velocity can also 
be associated with each region to model inter- 
granular variation. A suitable granular model 
is shown in table ^ Using this granular model, 
calculated spectral lines closely match observed 
spectral lines in width and asymmetry. 



Parameter 


Upflow 


Transition 


Downflow 


Area 


0.45 


0.40 


0.15 


Brightness 


1.07 


0.97 


0.87 


Weighting 


0.48 


0.39 


0.13 


V 


0.577 





-1.072 


v s 


368 


368 


368 


micro 


1.58 


3.67 


3.67 


macro 


1.6 


1.6 


3.5 



Table 1 



Granular model. All velocities are in km/s, and 
heights in km 



7. Spectral synthesis 

If we follow the steps outlined here, we can cal- 
culate the emergent spectrum of the sun. The 
matching of observed and calculated line pro- 
files allows us to investigate the processes respon- 
sible for the strength and shape of the line — 
abundance, damping, granulation, etc. Spectral 
synthesis, a valuable tool for studying the atmo- 
spheres of the sun and other stars, can be readily 
implemented on a typical PC. A number of spec- 
tral synthesis codes are available |12I13| . 
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